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regained significant interest with detailed spectral and timing studies of the X-ray 
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emission from microquasars, the recent detection by the HESS collaboration of 
very-high-energy 7-rays from the microquasar LS 5039, and the earlier suggestion 



O 

of jet models for ultraluminous X-ray sources observed in many nearby galaxies. 



Here we study the synchrotron and Compton signatures of time-dependent elec- 
tron injection and acceleration, adiabatic and radiative cooling, and different jet 
geometries in the jets of Galactic microquasars. Synchrotron, synchrotron-self- 
Compton, and external-Compton radiation processes with soft photons provided 
by the companion star and the accretion disk are treated. An analytical solution 
is presented to the electron kinetic equation for general power-law geometries of 
the jets for Compton scattering in the Thomson regime. We pay particular atten- 
tion to predictions concerning the rapid flux and spectral variability signatures 
expected in a variety of scenarios, making specific predictions concerning possi- 
ble spectral hysteresis, similar to what has been observed in several TeV blazars. 
Such predictions should be testable with dedicated monitoring observations of 
Galactic microquasars and ultraluminous X-ray sources using Chandra and/or 
XMM-Newton. 



Subject headings: gamma-rays: theory — radiation mechanisms: non-thermal - 
X-rays: binaries 
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1. Introduction 

The high-energy emission from X-ray binaries (XRBs) is generally believed to be pow- 
ered by the accretion of matter from a stellar companion onto a compact object of typically 
a few solar masses. In many sources, this accretion process is associated with the expulsion 
of collimated, mildly relativistic bipolar outflows (jets) most likely perpendicular to the ac- 
cretion disk. In a few sources, e.g., GRS 1915+105 (Mirabel & Rodriguez 1994), VLBI radio 
images have shown a spatially-resolved radio jet during episodes of quasi-steady and hard 
X-ray emission (Dhawan et al. 2000). 

In the standard picture, the high-energy (X-ray - 7-ray) spectra of X-ray binaries gen- 
erally consist of two major components: A soft disk blackbody with a typical temperature of 
kT ~ 1 keV, and a power-law at higher energies. Neutron-star and black-hole X-ray binaries 
exhibit at least two main classes of spectral states, generally referred to as the high/soft state, 
and the low/hard state (for a review see, e.g. Liang 1998; McClintock & Remillard 2004). 
The high-energy spectra of X-ray binaries in the soft state are characterized by a thermal 
blackbody component, believed to be associated with thermal emission from an optically 
thick, geometrically thin accretion disk (Shakura & Sunyaev 1973), and a power-law tail 
with a photon index T > 2. Generally, no high-energy cutoff of the high-energy power- law is 
detected. In the hard state, the spectrum is dominated by a power-law, with a slope T < 2 
and a cut-off at ~ a few hundred keV. 

While, in the conventional view of X-ray binaries (including microquasars), the X- 
ray and 7-ray emission is attributed to Comptonized emission (Sunyaev & Titarchuk 1980; 
Titarchuk 1994) arising from hot thermal (kT 3> IkeV) or relativistic, non-thermal electrons 
close to the black hole (e.g. Liang & Price 1977; Bisnovatyi-Kogan & Blinnikov 1977; Shapiro, 
Lightman, & Eardley 1976; Narayan & Yi 1994; Chen et al. 1995), the tentative EGRET 
detections of at least two Galactic microquasars at MeV -GeV 7-ray energies, namely LS 
5039 (Paredes et al. 2000) and LSI+61°303 (Gregory & Taylor 1978; Taylor et al. 1992; 
Kniffen et al. 1997), the detection of X-ray jet structures in several microquasars using 
Chandra and XMM-Newton (Corbel et al. 2002; Tomsick 2002, e.g.,), and, most recently, 
the detection of very-high-energy (VHE) 7-ray emission from LS 5039 (Aharonian et al. 
2005) have re-ignited interest in jet models for the high-energy emission from microquasars, 
analogous to the commonly favored models for blazars (for a recent review see, e.g. Bottcher 
2002). 

A jet origin of the X-ray emission of microquasars has previously been suggested by sev- 
eral authors, e.g., Markoff et al. (2001, 2003a,b), who discussed the possibility of synchrotron 
emission from relativistic electrons in the jet extending from the radio all the way into the 
X-ray regime. This idea is particularly well motivated by the tight temporal correlation 
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between the radio and hard X-ray emission in micro quasars, apparently anti-correlated with 
the soft X-ray emission (e.g. Corbel et al. 2000, 2001). 

Additional motivation for the investigation of jet models for the X-ray emission of mi- 
croquasars comes from the suggestion by Georganopoulos, Aharonian & Kirk (2002) that 
ultraluminous X-ray sources (ULXs, with X-ray luminosities L x > 10 39 erg s -1 ), detected 
in many nearby galaxies (see, e.g. Makishima et al. 2000; Fabbiano et al. 2001), may be 
microquasars viewed at very small angles with respect to the line of sight (hence dubbed 
"microblazars" ) . However, more recent observational evidence, including indications of pre- 
dominantly thermal spectra of ULXs (e.g. Miller et al. 2003) and the discovery of X-ray 
ionized nebulae around the optical counterparts of ULXs (Pakull & Mirioni 2003; Kaaret et 
al. 2004), provide increasing support for alternative models. 

Motivated by the new observational results summarized above, various authors have 
been working on the expected high-energy emission signatures from microquasars. In ad- 
dition to the jet-synchrotron models mentioned above, Bosch-Ramon et al. (2005a,b) and 
Romero et al. (2005) have investigated various scenarios for the leptonic or hadronic origin 
of high-energy emission from microquasars, aiming, in particular, at explaining the potential 
EGRET detections of LS 5039 and LSI+61°303. Those papers concentrate on a detailed 
description of the high-energy emission processes in a steady-state configuration. However, 
if leptonic emission mechanisms play a dominant role in the high-energy radiation of micro- 
quasars, then one would expect to see significant spectral variability, possibly similar to the 
spectral hysteresis observed in the X-ray spectra of several high-frequency peaked BL Lac 
objects, such as Mrk 421 (Takahashi et al. 1996) or PKS 2155-304 (Kataoka et al. 2000). 

Indications of such behavior have been seen in XTE J1550-564 (Rodriguez et al. 2003), 
but X-ray variability correlated with flaring behaviors detected at GeV energies with GLAST 
or at TeV energies with the ground-based air Cherenkov telescopes would provide strong 
evidence that leptonic processes are responsible for the 7-ray emission. A detailed time- 
dependent analysis of the high-energy emission signatures of jet models of microquasars has 
so far been restricted to a study of the X-ray timing signatures of an empirical power-law 
model for the high-energy emission (Kording & Falcke 2004), and a study of the expected 
QPO features as well as broad-band spectral variability correlations on longer time scales 
(Yuan et al. 2005). Reig, Kylafis, & Giannios (2003) (see also Giannios, Kylafis & Psaltis 
2004; Giannios 2005) proposed a model of X-ray variability from microquasar jets, in which 
the variability is driven by a variable soft photon source (the accretion disk), with a steady- 
state jet configuration. Those authors are focusing on predictions concerning the power 
density spectra as well as time and phase lag features, similar to previous work for accretion- 
disk corona models by Kazanas, Hua & Titarchuk (1997); Hua, Kazanas & Titarchuk (1997); 
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Hua, Kazanas & Cui (1998); Bottcher & Liang (1998) and Bottcher (2001). No predictions 
concerning X-ray spectral hysteresis in HID diagram form are made. 

Since it is likely that many different radiation components (synchrotron, synchrotron- 
self-Compton, external-Compton) are contributing to the high-energy emission from micro- 
quasars in the low-hard state, the rapid variability patterns expected in realistic jet models 
of microquasars might be more complex than the features previously investigated. In this 
paper, we present a detailed study of various plausible scenarios of electron injection and 
acceleration into a relativistically moving emission region in a microquasar jet, and subse- 
quent adiabatic and radiative cooling. We pay particular attention to the X-ray spectral 
variability, as motivated above. In §2, we present a general outline of the assumed model 
geometry and choice of parameters for a baseline model used for our investigation. This 
section also contains a brief discussion of the emission mechanisms which we include in our 
calculations. In §3, we present an analytic solution to the electron kinetic equation in the 
case of power-law jet geometries, which is used for our time-dependent study of the radiation 
signatures from our model. Starting from our baseline model, we then perform a parameter 
study to investigate the imprint of various parameter choices on the expected spectral and 
variability patterns from microquasars, which is presented in §4. We summarize in §5. 



The geometry of our model setup is illustrated in Fig. 1. The accretion flow onto the 
central compact object is ejecting a twin pair of jets, directed at an angle 9 with respect 
to the line of sight. Two intrinsically identical disturbances, containing non-thermal plasma 
(blobs) originate from the central source at the same time, traveling in opposite directions 
along the jet at a constant speed Vj = f3j c. Let d be the distance to the source, and 
fi = cos 9. The time at which any radio component is observed, is denoted by t Q bs- If we 
wish to distinguish between the approaching and the receding blobs, we do so using the 
subscript a for the approaching and r for the receding component of the jet. The intrinsic 
times (in the stationary frame of the central source) at which those components are being 
observed, are 



corresponding to linear distances of x a r = Vjt* &T from the central source. Here and through- 
out the remainder of the paper, times t* measured in the rest frame of the central source are 
denoted by a superscript '*', while times t without superscript refer to the co-moving frame 
of the blob. We assume that over a limited range in distance x Q < x < x 1 , relativistic elec- 
trons are accelerated and injected in the emission region. Most scenarios for the acceleration 



2. Model Setup 
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and injection of relativistic electrons in the jet will lead to a power-law spectrum of injected 
electrons, which we parametrize in terms of an electron injection function, 

- -r^ir^ = Q° "<i q ^ min ' ^ max ) H ^ to ' *0> ( 2 ) 

where H(y;y ,yi) is defined as 1 for y < y < y\ and elsewhere, and to,i = The 
normalization constant Qo for the electron energy distribution is related to the injection 
luminosity L in j through 

ifg^2; 
if q = 2. 



Q = <( m ^ 2 i^-iL q n ) ' ' ' (3) 



The blob's (transverse) radius, R±, scales with distance from the central engine as 
R± = R° ± (x/x ) a , i.e., a = corresponds to perfect collimation, and a = 1 describes a 
conical jet. In the following, we will consider values of < a < 1, and for the purpose of 
an analytical estimate, we assume no significant expansion along the direction of motion. 
Following the arguments given in Atoyan & Aharonian (1997), we choose a magnetic-field 
dependence on distance from the central black hole as B(x) = B (R ± /R°_ L )~ 2 = B (x/x )~ 2a . 
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Each electron injected into the emission region at relativistic energies, will be subject 
to adiabatic and radiative cooling, described by 



losses. The second term describes synchrotron and Compton losses, with u = UB+u ra a, where 
ub = B 2 /8n is the magnetic-field energy density, and u ra d is the seed photon energy density 
for Compton scattering in the Thomson regime. The term -u ra d consists of contributions from 
the X-ray emission from an optically thick accretion disk, from the intrinsic synchrotron 
radiation, and from external photons from the companion star. From the milli-arcsecond 
resolution of the VLBA monitoring observations which have resolved superluminal-motion 
components in Galactic microquasars, we can estimate that such components appear on 
linear scales of x ~ 1CT 3 x_ 3 pc with x_ 3 ~ 1 - 100. The energy density of disk photons 
in a point-source approximation is then u cxt ~ 3 x 10~ 5 L 38 (x_ 3 )~ 2 Tj 2 ergs cm~ 3 . The 
magnetic-field energy density of u B ~ 4 x 10~ 2 Bq ergs cm~ 3 where B G is the magnetic field 
in Gauss. Given the estimated magnetic-field values of ~ a few tenths of Gauss estimated 
on scales of several mpc for GRS 1915+105 (Atoyan & Aharonian 1999), the magnetic-field 
decay mentioned above suggests initial magnetic field values of -Bo ~ several thousand Gauss 
at injection sites x of a few thousand R g from the black hole. Since there does not seem to be 
any evidence for a "Compton catastrophe" in Galactic jet sources, we may assume that the 
soft photon energy density of the synchrotron radiation field is significantly lower than the 
magnetic-field energy density. However, in high-mass X-ray binaries (HMXB) like Cyg X-l, 
LS 5039, or LS I +61°303, it is possible that the energy density in the soft photon field from 
the companion star is dominant over the magnetic-field energy density. The soft-photon 
energy density from a companion star of luminosity = 10 39 L 39 ergs s _1 at a separation 
of s = 10 12 s 12 cm from the black-hole primary is u* ~ 3 x 10 3 (T 2 L^/y 2 2 ) ergs cm" 3 , where 
Vi2 = V x i2 + s i2 1S ^e distance between the blob and the star in units of 10 12 cm. 

The standard parameter choices for our "baseline" model, are broadly representative 
of GRS 1915+105 in the low/hard state, which also give equipartition between the energy 
densities of the relativistic electrons and the magnetic field in the ejecta. Repeated VLA 
observations of ejection events from radio emitting clouds of GRS 1915+105 have established 
that the ejecta move with a true speed of (3 — 0.92, giving a bulk Lorentz factor of Tj = 2.5, 
at an angle = 70° with respect to the line of sight (Mirabel & Rodriguez 1994). The recent 
measurement of rotational broadening of the early-type K giant companion star, combined 
with the orbital parameters lead to a black hole mass of 14 ± 4M (Harlaftis & Greiner 
2004). For our parameter study, we have adopted a black hole mass of 15M©. Radial 
systemic velocity and radio data calculations (Fender et al. 1999; Greiner et al. 2001) have 



d>y 1 dV B 7 4 




(4) 



dt V B dt 3 3 



in the co-moving frame of the blob, where the first term on the r.h.s. describes the adiabatic 
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yielded a distance of d ~ 12kpc to the source. Our reference values for the orbital separation 
of the binary and the donor star's luminosity follow from a mass function of 9.5 ± 3M Q 
derived by Greiner et al. (2001) from infrared observations. For the disk, a mass accretion 
rate of 1% in Eddington units, translating to ~ a few times 10 _9 M Q /yr for a 15M Q black 
hole, is a reasonable estimate in the low-hard state (McClintock & Remillard 2004; Belloni, 
Migliari & Fender 2000). We inject electrons into the jet over a distance of about 100-R 9 , 
and follow the evolution of the jet for about 5 days. The magnetic field values in the jet 
are hard to constrain, but we start at an initial value of B = 5000G at the base of the jet, 
resulting in B ~ 0.2 — 0.3G after ~ 4.5 days from the start of ejection (Atoyan & Aharonian 
1999). The above parameter choices are summarized in Table 1. 

For the purpose of an analytic treatment of the electron kinetic equation, Eq. 4, we 
re-write it in the form 



g?7 7 



-4a 



-2 



+ /.(*) 7 2 



(5) 



with m = 2a/3, 



- 4 CJj B % 
3 8ir m P c 2 ' 



(6) 

and f*(t) is a function describing the decline of the stellar radiation field energy density 
with distance from the black hole. For the purpose of an analytical description, we will 
approximate f*(t) as a broken power-law 



/.(*) = 



x 2 + s 2 



£ s 2 for x < s, 
£ x~ 2 for x > s 



(8) 



where £ = a T L^ T 2 /3irm e c 2 . An analytic solution to Eq. 5 will be derived in the 
Appendix and discussed in §3. 

Once a solution 7i(^;7, t) to Eq. 5 is found, the electron distribution at any point in 
time (and thus at any point along the jet) can be calculated through the expression 



N, 



t 

(r,t) = J 



dti 



d"fi dti 



t 



dji, 
d'y 



(9) 



Additional restrictions on the integration in Eq. 9 and possible solutions of 7(i;7i,^) are 
imposed by the Heaviside functions in Eq. 2 and will be specified in §3. 
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Radiation mechanisms included in our simulations are synchrotron emission, Comp- 
ton upscattering of synchrotron photons, namely synchrotron self-Compton (SSC) emission, 
and Compton upscattering of external photons. With the time-dependent (and thus in- 
dependent) non-thermal electron spectra in Eq. 9, we then use a 5-function approxima- 
tion to estimate the vF v synchrotron spectral output fl y at a dimensionless photon energy 
e = hu/m e c 2 (in the observer's frame): 

ft = -r^ (10) 



4-7T d\ m e c 2 dt de 



and 



d 2 E 4 



dt de 3 



-CCTtUbYN^) (11) 



f? (e, t obs ) = -D 4 i^f) 1% N e ( lsy , t) , 7sy ee J^- (12) 



giving 

fSyUf LA = F> ^'6nd 2 ) 'sy^ev/sy,v, <*y-}/ DeB 

(Dermer, Sturner, & Schlickeiser 1997), where D = [Tj (1 — f3jij)]~ l is the Doppler boosting 
factor, e B = B/B cv with _B cr = m 2 c 3 /he = 4.414 x 10 13 G, defined as the field at which the 
cyclotron quantum number equals the rest mass energy of the electron, t b s is the observer 
time, so that t obs = t/D = t* ai /DTj, and E is the total synchrotron energy radiated in the 
system. 

For the time-dependent vF v spectral output ff° due to Compton upscattering of ex- 
ternal photons from the star, we use the Thomson approximation and represent the star as 
a monochromatic point source. We find 



ff% ^^^) l!Wh ^^ (13) 

where u ext is the Doppler-boosted seed photon energy density from the companion star (in 
the co- moving frame), given by 

"- = D -4^V < 14 » 
where e ext is the dimensionless photon energy due to the companion star in the co-moving 
frame. It is related to the characteristic photon energy from the star e* = 2.7 /cT^/m e c 2 
through e ext = D*e*. T* is the effective surface temperature of the star, and = Tj(l — 
PrfJ'txt) Lorentz transforms e* into the comoving frame. 

Here we have omitted an angle-dependent factor that arises in a more accurate derivation 
for the spectrum of Thomson-scattered stellar radiation, which is important when the orbital 
variability of the star is considered (Dermer & Bottcher 2006). 
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The photon spectrum due to Compton upscattering of photons from the accretion disk 
is obtained by adopting a similar formalism as that used for the companion star. We approx- 
imate the disk as a thin annulus at the radius -R max where the differential energy output is 
maximized. We use the optically thick, geometrically thin, gas-pressure dominated accretion 
disk solution of Shakura & Sunyaev (1973) predicting a blackbody spectrum according to a 
temperature distribution TD(R max ) given by 

kT 

®(Rmax) — ~(Rmax) (15) 

m e c z 



3 



\M B ) \M e /yr) \ R„ J \ V ) 

where R g is the gravitational radius, -R max — 9.375 R g , Td is the corresponding disk temper- 
ature, M is the black hole mass in units of M , and M is the mass accretion rate. 

We also include SSC emission, keeping in mind that our analytical approach outlined in 
the following section will not be applicable to situations in which the radiative output (and 
thus electron cooling) from SSC dominates the bolometric luminosity, since we are neglecting 
SSC cooling of the electrons. We use the formalism of Tavecchio et al. (1998) to calculate 
the SSC spectral output in the Thomson regime. In order to be able to use this formalism, 
we approximate our time-dependent population of relativistic electrons at any given time by 
a broken power law distribution with indices n\ and n 2 and break Lorentz factor 7^: 

N{l) = l Kll ' ni if7<76 (17) 

where the coefficients K\ and K 2 are related through the kinetic condition K\^ ni = K 2 % n2 . 
The value of 7^ is determined by the solution of the kinetic equation for electron cooling. 
For weakly beamed jet emission, some fraction of the radiation is expected to impinge upon 
the cooler material in the accretion disk and lead to fluorescent Fe Ka line emission and a 
reflection hump (Beloborodov 1999). The amount of radiation irradiating the disk depends 
on the distance between the jet emission zone and the disk, the jet inclination angle, and 
the bulk Lorentz factor Tj. X-ray emission originating at distances ~ 10 3 R g from the disk, 
together with (3j = 0.92 in our model, corresponds to the case of a "synchrotron dominated 
jet" (Markoff & Nowak 2004; Markoff, Nowak & Wilms 2005), where the fraction of re- 
flected emission has been calculated to be only l%-2%. We therefore do not include this 
component in our present calculations. 
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3. Analytic Solutions to the Electron Kinetic Equation 

The dynamic equation (5) can be solved analytically for arbitrary values of a, as long 
as the distance dependence of the external radiation fields, characterized by /*(£), can be 
piecewise approximated as a power-law. The solution to the electron cooling problem, as 
derived in the Appendix, is given by 

7 (t; 7i , U) = t~ m \— + ^ Ul-m-Aa _ t l-m-4a\ + 

d ' t - {1+m) - t-< 1+m >) + f u{t') f- m df] 1 (is) 



l + m 

where (5 = u Q t^ a and 5 = v-st^. Note that one can easily incorporate as many additional 
external photon components with different decay slopes and normalizations as required for 
any specific problem at hand. 

The evaluation of the time-dependent electron distribution becomes more easily 
tractable if we transform the ti integration in Eq. 9 into an integral over the injection energy 

10 



7i 



i .max 



N e (l,t) = Q J dlill 



d 7i 



dU 



(19) 



7i,min 

The expression for the Jacobian in Eq. 19 based on our analytical solution is given in the 
Appendix (Eq. A4). The boundaries of the integral (19) follow from the Heaviside functions 
in Eq. 9, which also provide the relevant limits on the values of 7 for which the electron 
distribution is non-zero. Depending on whether the emission region is currently within the 
zone of electron injection (t < ti) or beyond it (t > ti), we find two sets of boundary 
conditions. 

In the case t < t±, we have 

7 i,min = 7 , ^Q) 

7 i,max = min ( 7max , 7i [i ; 7 , t}) , 
with non-zero electron distribution for 

l(t\ 7min , t ) < 7 < 7max (21) 
where 7i [to; 7 ,t] is found by inversion of solution in Eq. 18. 
In the case t > t±, we find 

7ijmin = max ( 7min , 7i [i J ; 7 , t] ) , 
7 i,max = min ( 7max , 7i [t ; 7> *]) » 
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with non-zero electron distribution for 

7(*l 7min, t ) < 7 < 7(*l 7max, (23) 

The solution, Eq. 18, reduces to a particularly simple form in the limit of negligi- 
ble radiative losses, in which the electron kinetic equation reduces to 7 = — 771(7 /*)> an d 
7(t;7i,ti) = 7; (t/ti)~ m . The local electron spectrum is then given by 

t- \^ ft- ^ v 



N: di (r,t) = Q,i- q t^^ (24) 

77 

where ti jm i n and ti imax are determined by the Heaviside function in Eq. 2, and 77 = (q— 1) m+1. 
In the limit t 3> ti and 7 <C 7 ma x(i), the local electron spectrum reduces to 

JV^( 7 ;t)«^7-« (25) 



4. Results 

A large number of simulations have been performed to study the effects of the various 
model parameters on the resulting broadband spectra, light curves, and X-ray hardness 
intensity diagrams (HIDs). We start our parameter study with a baseline model for which we 
have used the standard model parameters discussed in § 2 and listed in Table 1. Subsequently, 
we investigate the departure from this standard set-up by varying (1) the initial magnetic 
field Bo, (2) the luminosity of the companion star L*, (3) the injection electron spectral index 
q, (4) the low-energy cutoff 7l of the electron injection spectrum, (5) the high-energy cutoff 
72 of the electron injection spectrum, (6) the injection luminosity L in j, (7) and the observing 
angle 9 obs and thus the Doppler boosting factor. The parameters used for the individual runs 
are quoted in Table 2. 

In Figure 2, we have compiled for our baseline model (a) a sequence of comoving electron 
spectra, (b) snapshot SEDs, and the time averaged photon spectrum, (c) light curves at 
various photon energies, and (d) a snap-shot spectrum, indicating the individual emission 
components. Fig. 2a illustrates the gradual build-up of the electron density in the emission 
region and their subsequent radiative and adiabatic cooling. Fig. 3 illustrates the electron 
cooling time scales (a) as a function of electron energy at the end of the injection period, and 
(b) as a function of time for a fixed electron energy of 7 = 10 3 . At the end of the injection 
period, the radiative cooling timescale is shorter than the elapsed time for electron energies 
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of 7 > 7b = 2.4 x 10 3 . In our baseline model, the magnetic field energy density is higher 
than that of the external photon field, with 

^1 „ 2D * L * (26) 
u B B 2 y 2 c V } 

attaining a value of 0.36 at the end of the injection period. Because of the rapid decline of 
the magnetic field, B oc x~ 2a oc x~ 0M , the cooling time scale increases with time in our base 
model, as t sy oc t 4a oc t 12 . For this reason, one can see a synchrotron cooling break in the 
electron spectra around % ~ 10 3 at early times, while, at the end of the injection period, 
this break blends in with a gradual high-energy cut-off and is no longer discernible. After 
the end of the injection period, synchrotron cooling is initially still dominant, but adiabatic 
cooling is gradually taking over as the dominant energy loss mechanism (see Fig. 3b). 

The time-dependent photon spectra in Fig. 2b indicate that the synchrotron and the 
Compton (disk) components decay even within the injection interval due to decaying source 
fields. The disk photon energy density begins to decline (oc t~ 2 ) immediately after the onset 
of injection, while the stellar photon field remains basically constant for x < s. This manifests 
itself in the gradual dominance of the EC (star) component over the EC (disk) component in 
its time evolution. The EC (disk) component produces the highest-energy photons, extending 
out to energies beyond ~ 100 GeV, but only during a very short, sub-second flash very 
close to xq. If the electron injection does actually start as close to the black hole as x ~ 
lOOORg , as we have assumed here, such high-energy radiation might suffer substantial 77 
absorption in the intense radiation field of the stellar companion (Moskalenko et al. 1993; 
Moskalenko & Karakula 1994; Bednarek 2000; Bottcher & Dermer 2005) and might thus not 
be directly observable. However, as noted in the introduction, Compton scattering of the 
stellar radiation field or possibly also hadronic processes might very well produce VHE 7-rays 
at larger distances from the central compact object (Romero et al. 2003; Bosch-Ramon et 
al. 2005a; Dermer & Bottcher 2006) which could escape the compact radiation field of the 
stellar companion without significant 77 attenuation. The radio spectrum is usually flat to 
inverted and extends to the infrared and even beyond (Fender et al. 2000). The spectral 
cut-off above the GHz range in Fig. 2b is a result of our choice of 71 = 10, the lower energy 
cut-off of the electron distribution. Effects of varying 71 are investigated in Section 4.5, 
showing a spectral shift towards lower frequencies with lower values of 71. The short-time 
X-ray spectral variability predictions, which are the main focus of this paper, are only weakly 
dependent on the choice of 71, at least within a realistic range of 71 ~ 1 — 10. 

We also note (see Fig. 2d) that the contribution from the SSC photons to the total 
photon spectra is negligible in our base model. 

Figure 2c shows light curves at different photon frequencies, one in the UV, two in X- 
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rays, and two at 7-ray energies. From the light curves and the snap-shot spectra in Fig. 2b, 
one can see that the synchrotron-dominated optical - UV spectra reach their maximum very 
rapidly and remain at an approximately constant level until the higher-energy light curves 
reach their maxima, around the end of the injection period. The approximately constant 
level of the optical - UV radiation results from the opposing effects of an increasing electron 
density and a decreasing magnetic field. There is no significant time delay between the light 
curve peaks at the various X-ray and 7-ray energies. At 24 MeV, the dip in the light curve 
before the end of injection is a result of the very rapid decay of the photon field density from 
the disk. Here, the EC (disk) component dominates very early on, but declines rapidly and 
becomes dominated by the rising EC (star) component within less than 1 s. In general, a 
rapid change in temporal slope, as seen in various light curves in Fig. 2c, indicates a transition 
from one radiation component to another, passing through the observing range. It should be 
mentioned here that the hard low-frequency cut-off apparent in our photon spectra instead 
of a vFl yn oc z/ 4 / 3 behavior at v < ^(71) is a result of the 5-function approximation in our 
calculations. 

In the following parameter study, we focus on the time-averaged photon spectra, light 
curves, and X-ray HIDs, and explore the effect of variations of individual parameters on 
these aspects. Throughout our study, the dominant radiative cooling mechanism will either 
be synchrotron or EC (star). We find that for plausible choices of parameters the contribution 
from SSC cooling will always be negligible, as required by our analytical approach to the 
solution of the electron kinetic equation given in Eq. (18). 

4.1. Initial Magnetic field 

In this subsection, we investigate how different choices of the initial magnetic field Bo 
(at x ) influence the shape of the time-averaged photon spectra, light curves, and X-ray 
HIDs. The results are illustrated in Fig. 4. 

With increasing values of B , obviously, the total output in synchrotron emission in- 
creases, at the expense of the Compton components. As all other parameters remain con- 
stant, this also leads to a more rapid radiative cooling of the electrons, resulting in a cooling 
break at lower electron energies and, consequently, lower vF v peak frequencies of all radia- 
tion components. Notice that the time-averaged flux at ~ 10 keV remains almost unchanged 
for different magnetic field values. For the parameters adopted here, the X-ray emission is 
generally dominated by the EC (star) component, except for very high magnetic field values, 
in which case the synchrotron component initially extends into the X-ray regime. This is 
reflected in the very high X-ray fluxes in the highest-magnetic- field cases (see Fig. 4b). Note 
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that in the Compton-dominated X-ray regime (3.5 keV), the effect of different magnetic- 
field values on the light curve is only very minor. This is because this part of the spectrum 
is dominated by the low-energy end of the EC (star) component, produced by low-energy 
electrons, whose cooling is dominated by adiabatic cooling (see Fig. 3a), irrespective of the 
magnetic field. 

The shorter radiative cooling time scale with increasing B Q (and thus decreasing < 7 2 >) 
also leads to an even further decreasing contribution of the SSC component since 

~< 7 2 > a T n e R. (27) 

J-'syn 

Panels c) and d) compile the hardness intensity diagrams for the different values of 
Bq. Here and throughout the rest of the paper, we define the hardness ratio as the ratio 
of the X-ray fluxes at 2 - 10 keV (hard X-ray flux) to 0.1 - 2 keV (soft X-ray flux). The 
abrupt breaks in some of these tracks are an artifact of the 5 approximations used in some 
of our spectral calculations. The HIDs trace out characteristic hysteresis loops, changing 
their orientation. This behaviour is expected if the main contribution of the flux in a given 
energy bin transits between the high-energy end of one emission component to the low- 
energy end of another one (in this case, between the synchrotron and EC (star) components; 
see Bottcher & Chiang 2002). At X-ray energies around the synchrotron cut-off, the flux 
maxima occur at significantly different hardness ratios for different values of B (Fig. 4d). 
Specifically, within the synchrotron-dominated, clockwise spectral hysteresis loops, the peak 
flux is accompanied by a harder spectrum for higher magnetic fields. This is a consequence 
of the increasing frequency of the synchrotron peak with increasing magnetic field. Within 
the Compton-dominated, counterclockwise spectral hysterersis loops, the peak hard X-ray 
flux is accompanied by a softer hardness ratio for larger values of B , which is a consequence 
of the more efficient electron cooling. 



4.2. Luminosity of the Companion Star 

In our base model, as in most cases with substantial mass transfer rates onto the compact 
object, the external source photon field for Compton scattering is dominated by the star 
light of the stellar companion. For this reason, we performed a series of simulations with 
increasing values of the stellar luminosity to investigate the influence of a varying external 
photon density on the broadband spectra, light curves, and spectral variability patterns. 
The results of these simulations are illustrated in Fig. 5. 

The time averaged photon spectra are shown in Fig. 5a, which clearly shows the emer- 
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gence of the EC (star) component at X-ray and 7-ray energies. The impact of an increasing 
star luminosity on other spectral components (at radio, optical, UV, and GeV 7-rays) is 
small as long as the synchrotron cooling time-scale is shorter than the EC (star) cooling 
time-scale (solid and dotted curves). When EC cooling becomes dominant over synchrotron 
cooling, it leads to a reduced power output in the synchrotron and EC (disk) components as 
well as a shift of the vF v peak frequency of all radiation components towards lower values 
(dashed and dot-dashed curves). 

Fig. 5b shows that the impact of a strong EC (star) component is quite prominent 
in the light curve at keV energies. Along with an increase in the peak vF v flux (Compton 
peak), an increase in the star luminosity also leads to a faster decay of the light curves in 
the high-luminosity cases, where the star light radiation energy density dominates over the 
magnetic-field energy density. 

The HIDs are plotted in Fig. 5c, d and indicate X-ray spectral hysteresis. For low 
values of L*, the X-ray fluxes are initially dominated by the synchrotron flux, leading to 
very soft spectra at the time of the peak soft flux. As the flux decays after the end of 
the injection period, the EC (star) component begins to sweep through the X-ray bands, 
producing increasingly hard spectra during the decaying portions of the light curves (i.e., 
the low-flux portions of the HIDs). As L* increases towards very high values, a secondary, 
counterclockwise loop in the soft-flux HID (Fig. 5c) emerges when the EC (star) component 
sweeps through the X-ray bands. In that case, the EC (star) component dominates the 
hard X-ray flux very early on, causing the overall pattern of the hard flux HIDs (Fig. 5d) 
to consist only of this counterclockwise loop. We conclude, that the observation of such a 
counterclockwise spectral hysteresis is a diagnostic of the dominance of the EC component 
in the X-ray regime. 



4.3. Electron Injection Luminosity 

The effect of an increasing injection luminosity, corresponding to a higher density of 
injected relativistic particles in the emitting region, is illustrated in Fig. 6. This leads to a 
corresponding increase in the overall bolometric luminosity and a stronger energy output in 
the SSC component, relative to the synchrotron, as given in Equation 27. 

As long as the SSC emission is not dominant, this essentially only leads to a higher 
flux level in the spectra and light curves. At high values of L in j, the SSC emission begins to 
play a non-negligible role in the X-ray regime, which leads to a slightly shallower slope of 
the X-ray light curves during the rising phase (t < ti). Also, for distances x < s, the decay 
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slope of the X-ray light curve is steeper for higher values of L in j . Both of these effects are a 
consequence of the fact that the magnetic field decay is faster than that of the star photon 
field, and as a consequence, the decay in seed photon density for the SSC emission is faster 
than that for the EC (star) emission. 

Panels c and d of Fig. 6 show the various HID tracks for different injection luminosities. 
The main effect at the soft band HIDs (Fig. 6c) is visible in a slight hardening of the 
spectra during the decay phase for high values of L in j as the SSC emission gradually begins 
to play a non- negligible role. In the hard band HIDs (Fig. 6d), one can see that for low 
Li n j, the luminosity peak is reached very early on, and is dominated by the synchrotron 
component leaking into the 2-10 keV band. A secondary counterclockwise loop emerges 
as the EC (disk) component sweeps through the 2-10 keV band. As L in j increases (along 
with the SSC contribution), the secondary loop extends towards higher relative flux values, 
until, for the highest value of L in j, the peak flux occurs as the HID tracks goes through the 
counterclockwise hysteresis loop. There is no significant difference in the local spectral index 
(or the hardness ratio) at the time of peak flux of the Compton-dominated hysteresis loop. 



4.4. Injection Electron Spectral Index 

The effect of changes in the electron injection spectral index q manifests itself obviously 
in the photon indices of the broadband spectra. Fig. 7a shows the time-averaged photon 
spectra for different values of q. We see that this spectral change results in a shift of the 
radiation peaks toward higher frequencies as the injection spectrum hardens. We see a 
decreasing value of the peak flux for higher values of q (softer injection spectra) because of 
a decrease in the density of higher energy electrons with increasing injection spectral index. 
This causes a larger fraction of the injected energy to go into adiabatic rather than radiative 
losses (see Fig. 3a). Therefore, the jet becomes radiatively less efficient with increasing q. 
This also goes in tandem with a decreasing contribution from SSC, which follows directly 
from Eq. 27. 

As illustrated in the light curves in Fig. 7b, a higher value of q leads to a steeper decline 
in flux once electron injection stops. At the energy of 3.55 keV, the primary contribution 
is from the low-energy end of the Compton component, which is increasing with increasing 
spectral index. For this reason, as a higher fraction of the injected particle energy is in- 
jected near the low-energy cutoff of the electron distribution for increasing values of q, the 
peak medium-energy X-ray flux during the early phase of the simulated flares is positively 
correlated with q. During the later decay phase, the light curves converge to rather similar 
shapes and flux levels, though the decay slope is still positively correlated with the spectral 
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index, as expected. 

The X-ray spectral hysteresis at energies just below and above the synchrotron cut- 
off is shifted according to the change in electron injection spectral index, but the overall 
characteristics remain unchanged (see Fig. 7c). The peak flux is reached at a higher value 
of the hardness ratio for a softer injection spectrum. The hardness ratio decreases during 
the rising portion of the soft X-ray flux light curve, and increases during the decaying phase 
as the synchrotron component transits through the soft X-ray bands and gives way to an 
increasing contribution of the Compton components. The hard X-ray flux is initially also 
dominated by synchrotron emission, except in the case of the softest injection index, q = 2.8. 
As the synchrotron component rapidly moves out of the hard X-ray band, accompanied by 
a decreasing hardness ratio, the Compton component becomes dominant, and a secondary, 
counterclockwise hysteresis loop develops in the cases of soft electron injection spectra. 



4.5. Low-Energy Cutoff of the Electron Spectrum 

Since we have parameterized our injection spectra through an injection luminosity and 
the bulk of the kinetic energy of electrons in the jet is carried by the lowest-energy electrons 
for spectral indices of q > 2, a change of the low-energy cutoff of the electron spectrum, 7 1; 
naturally has quite dramatic effects on the radiative signatures from microquasar jets. These 
are illustrated in Fig. 8. Because of the dominance of adiabatic losses at low electron energies, 
the jet becomes radiatively more efficient for higher values of 71. Consequently, the total 
radiative energy (fluence) increases, as also seen in the previous section as a consequence of a 
harder electron injection spectrum. This also goes in tandem with an increased contribution 
of SSC to the high-energy emission. A higher value of 71 shifts the low-frequency cutoff 
of the individual radiation components towards higher frequencies. However, since the vF v 
peak frequency is dominated by emission from the highest-energy electrons for q < 3 (in the 
slow-cooling regime), a change of the value of 71 does not affect the vF v peak frequency of 
the individual radiation components. The shift of the high-energy peak into the 7-ray regime 
for 71 = 10 3 as seen in Fig. 8a is a consequence of the SSC component becoming dominant 
over the EC (star) component compared to the lower values of 71. 

The increasing contribution from SSC for 7x = 10 2 and 7x = 10 3 is also reflected in the 
light curves plotted in Fig. 8b. It is evidenced in the faster onset, but then flatter rising 
portion of the flux for 71 = 10 2 and the faster initial decay of that light curve, untill the EC 
(star) component begins to move through the observing frequency (3.55 keV), resulting in a 
similar decay as in the lower-71 cases. In the case of 71 = 10 3 , the SSC component always 
dominates the 3.55 keV X-ray light curve. 
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The HIDs shown in Fig. 8c, d confirm these findings: In the synchrotron-dominated 
portions of the X-ray flux (soft X-rays), the spectral hysteresis loops are clockwise and 
rather similar; however, for higher values of 7 1; the hardness ratio becomes significantly 
lower because of the increasing low-energy cutoff of the Compton spectra. The HIDs for the 
hard X-ray flux show a similar trend as described in the previous section: In particular, for 
low or moderate values of 71, a spectral softening during the early flux decay is followed by a 
counterclockwise secondary hysteresis loop due to the influence of the EC (star) component. 

4.6. High-Energy Cutoff of the Electron Spectrum 

A higher cutoff of the electron injection spectrum will manifest itself in the broadband 
spectra primarily by the extension of all radiation components towards higher energies, as can 
be seen in Fig. 9a. Since only a minor portion of the overall particle energy is carried by the 
highest-energy particles, the overall energetics of the particle distribution and the bolometric 
luminosity and radiative energy output of the emission region and the individual radiation 
components remain virtually unchanged. For high values of 72, the individual radiation 
components are increasingly overlapping, leading to a much smoother overall shape of the 
broadband spectrum. 

The 3.55 keV light curves, plotted in Fig. 9b clearly show the marked difference between 
the case of relatively low 72, where the influence of the synchrotron component at 3.55 keV is 
very minor, and the high-72 cases, in which the synchrotron component extends significantly 
beyond that energy. For our standard choice of B Q = 5 kG and D = 0.583, the critical 
electron Lorentz factor for synchrotron emission at 3.55 keV is 7 syi o = 10 4 , while at the end 
of injection this value has increased to 7 sy; i = 4 x 10 4 due to the corresponding decline of 
the magnetic field over the injection length. The (co-moving) synchrotron cooling time scale 
for electrons emitting 3.55 keV synchrotron radiation at the end of the injection period is 
t sy = 0.2 s. Consequently, the synchrotron contribution to the 3.55 keV light curve disappears 
virtually instantaneously after the end of injection, so that the three light curves are basically 
identical beyond that point. 

Fig. 9c reflects the fact that during the rising phase of the light curve, the soft X-ray flux 
is clearly dominated by synchrotron emission. In the high-72 cases, the synchrotron cutoff 
is located beyond the soft X-ray regime, leading to a relatively hard soft-X-ray spectrum 
that gradually softens due to the development of a cooling break moving through the 0.1 
- 2 keV flux. For 72 = 10 4 , the effect of radiative cooling is immediately visible even in 
the soft X-ray regime, leading to a systematically softer hardness ratio near the time of soft 
X-ray peak flux. The secondary hard X-ray spectral hysteresis visible in the 72 = 10 4 case is 
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not evident in the higher-72 cases just because of the much higher peak flux, so that those 
hysteresis loops occur at much lower relative flux values (compared to the maximum flux). 

4.7. Doppler Boosting Factor 

We have investigated the effect of changing the Doppler factor D on the spectra, light 
curves and HIDs by changing the observing angle 6 obs . Fig. 10a illustrates the spectral 
change resulting from a varying observing angle, given essentially by u b s oc D and vF v oc D 4 
for intrinsically isotropic emission (note, however, that the dependence is somewhat stronger 
for EC emission, e.g., Dermer 1995). 

The 3.55 keV light curve (Fig. 10b) and HIDs (Fig. 10c,d) illustrate the gradual shift 
of the synchrotron emission out of the X-ray regime for lower Doppler factors (between 
individual runs), and a decaying magnetic field as well as the development of the cooling 
break, moving through the X-ray energy range (within individual runs). In particular, this 
leads to synchrotron-dominated clockwise spectral hysteresis loops in the soft X-ray HIDs, 
and the emergence of secondary, EC (star) dominated, counterclockwise hysteresis loops in 
the hard X-ray HIDs for lower beaming factors. At times when the synchrotron component 
does no longer have any influence on the X-ray light curves, the X-ray emission at 3.55 keV is 
dominated by the EC (star) component, and we observe two separate regimes, depending on 
the Doppler factor. For smaller Doppler boosting factors, the spectral break corresponding 
to star light scattered by electrons that were injected at 7 = 71 occurs always below 3.55 keV. 
In that case, the peak in the X-ray light curve corresponds to the end of electron injection 
at £b,obs — ti/D. Consequently, it shows up at earlier observed times for increasing values 
of D. For higher values of the beaming factor, initially the observing frequency is below the 
break frequency corresponding to 71. In that case, the X-ray light curve keeps rising even 
after the end of the injection period, and the peak corresponds to the time where the break 
frequency becomes lower than the observing frequency. Since this happens well after the end 
of electron injection, the time of peak flux in this case is in fact later than that for smaller 
beaming factors, which explains the later light curve peak for 9 = 5° compared to 9 = 20°. 
Since the cooling of electrons injected at 71 < 100 is dominated by adiabatic losses, their 
energy decreases as 7(2) oc t~ rn . Consequently, we expect a dependence tt>,obs oc _D _0 3//a once 
the transition to this regime has occurred. 
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5. Summary and conclusions 

We have presented a detailed parameter study of the spectra resulting from time- 
dependent injection and acceleration, and adiabatic and radiative cooling of nonthermal 
electrons in the jets of Galactic microquasars. Jet models of microquasars have recently 
attracted great interest, especially after the detection of VHE 7-ray emission from the high- 
mass X-ray binary and microquasar LS 5039, in combination with the tentative identification 
of several microquasars with unidentified EGRET sources. These detections have confirmed 
the idea that microquasars are a distinctive class of high and very high energy 7-ray sources. 

It remains an open question whether the high-energy emission from microquasars is 
associated with leptonic or hadronic primaries. Predictions for the correlated variability at 
X-ray and high-energy 7-ray energies provides a discriminant between the two possibilities. 
As shown here, X-ray hysteresis diagrams are predicted in nonthermal leptonic models of 
microquasars, similar to the situation for blazars. Detection of such variability behavior, par- 
ticularly if correlated with 7-ray flaring behavior as described here, would provide evidence 
in favor of leptonic models. 

Whereas previous studies have largely focused on spectral fits to the broadband emis- 
sion of microquasars in a steady-state approximation, we have focused on the spectral vari- 
ability features expected in generic microquasar jet models with electron injection and/or 
acceleration over a limited amount of time and length along the jet, representative of, e.g., 
internal-shock models. We have conducted a detailed parameter study, investigating the im- 
pact of variations of several fundamental model parameters on the broadband SEDs, X-ray 
light curves, and the rapid X-ray spectral hysteresis phenomena expected to arise in these 
scenarios. We have provided an analytical solution to the electron kinetic equation, taking 
into account radiation signatures of synchrotron, external Compton (with seed photons from 
the companion star and the accretion disk), as well as synchrotron self-Compton emission. 
In order to be able to work with an analytical solution to the electron kinetic equation, 
we restricted our analysis to Compton scattering in the Thomson regime. Consequently, 
our results are important to derive diagnostics of microquasar jet emission in the X-ray 
regime, to contrast predictions of thermal Comptonization models of these sources, but will 
be primarily applicable to microquasars that are not candidates of VHE 7-ray emission. We 
have neglected the angle dependence of the stellar radiation field (Bottcher & Dermer 2005; 
Dermer & Bottcher 2006), which would lead to a modulation of the 77 absorption trough 
as well as the Compton scattered stellar radiation spectra depending on the orbital phase 
(0o — showing the most absorption but highest Compton flux). The spectrum is expected 
to be further modified by pair cascades, redistributing some of the VHE radiative power 
to lower frequencies, and significantly increasing the transparency of the source (Aharonian 
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et al. 2005). A study incorporating full Klein-Nishina effects on Compton scattering, the 
angle dependence of the stellar radiation field and cascade processes will be presented in a 
forthcoming publication on this subject. 

Obviously, various spectral components (synchrotron, SSC, external-Compton) could be 
easily distinguished if detailed snapshot SEDs could be measured for microquasars, on the 
(often sub-second) time scales of their X-ray variability. Unfortunately, such detailed snap- 
shot broadband spectra are currently not available, and might not be available in the near 
future. Therefore, we have exposed several other features pertinent to the transition between 
different spectral components which will be more easily observable in realistic observational 
data of microquasars: 

• A sudden increase of a light curve slope at a fixed observing frequency, not accompanied 
by significant flaring activity at other wavelengths, usually indicates the passing of a 
new spectral component through the fixed observing frequency range. Most notably, 
this diagnostic can be used to investigate the presence of one or more external-Compton 
component (s) in the X-ray / soft 7-ray regime. 

• Clockwise spectral hysteresis in the hardness-intensity diagrams indicates the domi- 
nance of synchrotron emission (in particular, before the end of the injection period in 
our generic model setup). In this case, the frequency- dependent light curve decay will 
be a useful diagnostic of the magnetic field strength in the jet (e.g., Takahashi et al. 
1996). 

• Counterclockwise spectral hysteresis in the hardness-intensity diagrams indicates the 
dominance of Compton emission (similar to the case of blazars, see, e.g., Bottcher & 
Chiang 2002). 

• In our study, we found, quite often, a co-existence of clockwise and counterclockwise 
X-ray hysteresis loops, which would provide a particularly powerful diagnostic, as it 
would allow to probe the characteristic transition energy between synchrotron and 
Compton emission, and its time dependence. 

There are only very few papers which present X-ray variability in the form of hardness 
intensity diagrams. Those papers that do, are presenting long-term variability information 
from RXTE PCA observations, spanning time scales of hundreds of days, thus representing 
the spectral characteristics of state transitions rather than the short-term variability patterns 
that are the focus of our analysis (see, e.g. Fender, Belloni & Gallo 2004; Belloni et al. 2005). 
In such a representation, possible hysteresis on short (intraday) time scales are very hard 
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to identify. An example where such short-term spectral hysteresis may be identified, can be 
found in Homan et al. (2005: http:/ /tahti. mit.edu/opensource/1655/), which shows an HID 
for GRO 1655-40. Here, within the overall long-term variability pattern, individual short 
flares (over a few days) can be discerned. One such outburst around MJD 53600 shows a 
spectral hysteresis, which when transformed to the orientation of the axes of HIDs in this 
paper, results in a counter-clockwise loop, indicating Compton dominated emission from the 
outflow for this outburst. 

We conclude therefore that the X-ray variability as predicted by our model can be used 
as a diagnostic to gain insight into the nature of the high energy emission in microquasar jets. 
In particular, a transition between clockwise and counter-clockwise spectral hysteresis would 
allow not only the distinction between different emission components, but also parameters 
such as the magnetic field, the Doppler boosting factor, and the characteristic electron 
injection / acceleration time. 

The work of S.G. and M.B. was supported by NASA through XMM-Newton GO grant 
no. NNG04GI50G and INTEGRAL theory grant NNG05GK59G. The work of C.D.D. is 
supported by the Office of Naval Research and GLAST Science Investigation no. DPR-S- 



The general electron cooling equation (5) can be solved analytically by virtue of the 
substitution 
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A. Analytic Solution of the Electron Kinetic Equation 
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The Jacobian in Eq. (19) is 
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resulting in an electron spectrum at any given point along the jet: 
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Table 1. Parameter choices for our baseline model 



Parameter 


Symbol 


Value 


Black-hole mass 


M 


15 Mq 


Distance 


d 


3.75 x 10 22 cm 


Jet inclination angle 


0jet 


70° 


Bulk Lorcntz factor 




2.5 


Binary separation 


s* 


10 12 cm 


Luminosity of companion star: 




8 X 10 37 ergs s" 1 


Surface temperature of the companion star 


T, 


3 x 10 4 K 


Initial blob radius 


Ro 


10 3 Rg 


Jet collimation parameter 


a 


0.3 


Accretion Fraction 


m = M/M Edd 


0.01 


Accretion disk luminosity 


L D 


1.9 X 10 37 ergs s" 1 


Electron injection spectrum, low-energy cutoff 


Tmin 


10 


Electron injection spectrum, high-energy cutoff 


7max 


10 4 


Electron injection spectrum, spectral index 


1 


2.4 


Beginning of electron injection zone 


x 


10 3 Rg 


End of electron injection zone 


x l 


10 5 Rg 


Magnetic field at x 


Bo 


5 x 10 3 G 


Injection luminosity 




4.4 x 10- 5 L Edd 
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Table 2. Parameter choices for our parameter study 



Sequence 


B [G] 


L* [ergs s 1 ] 


1 


7i 


72 




#obs [°] 


D 


(1) 
(1) 
(1) 
(1) 
(1) 


1.0 x 10 3 
2.5 x 10 3 
5.0 x 10 3 
7.5 x 10 3 
1.0 x 10 4 


8 x 10 37 


2.4 


10 


10 4 


4.4 x 10~ 5 


70 


0.583 


(2) 
(2) 
(2) 
(2) 


5.0 x 10 3 


8 x 10 36 
8 x 10 37 
8 x 10 38 
8 x 10 39 


2.4 


10 


10 4 


4.4 x 10 -5 


70 


0.583 


(3) 
(3) 
(3) 
(3) 


5.0 x 10 3 


8 x 10 36 


1.2 
1.8 
2.4 
2.8 


10 


10 4 


4.4 x 10~ 5 


70 


0.583 


(4) 
(4) 
(4) 
(4) 


5.0 x 10 3 


8 x 10 36 


2.4 


1 

10 
100 
1000 


10 4 


4.4 x 10~ 5 


70 


0.583 


(5) 
(5) 
(5) 


5.0 x 10 3 


8 x 10 36 


2.4 


10 


10 4 
10 5 
10 6 


4.4 x 10~ 5 


70 


0.583 


(6) 
(6) 
(6) 
(6) 


5.0 x 10 3 


8 x 10 36 


2.4 


10 


10 4 


2.2 x 10 -5 
4.4 x 10" 5 
2.2 x 10~ 4 
4.4 x 10~ 4 


70 


0.583 


(7) 
(7) 
(7) 
(7) 


5.0 x 10 3 


8 x 10 36 


2.4 


10 


10 4 


1 x 10 35 


5 

20 
35 
70 


4.60 
2.88 
1.60 
0.583 
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Fig. 2. — The baseline model (for parameters see Table 1): (a) Time-dependent electron 
spectra, beginning at xo, increasing by a multiplicative step of 10 between individual curves. 
The arrows indicate the sense of time evolution; the last step corresponds to a co-moving 
time of 87.6 s. (b) Time-dependent vF v photon spectra, corresponding to the same times 
as shown in panel a. The heavy gray curve is the time-integrated flux (fluence, in units of 
Jy Hz s). (c) Light curves at five frequencies, as indicated in the legend, (d) Snap-shot 
spectrum just before the end of the injection period, corresponding to t ~ 1.5 s, showing the 
individual emission components. 
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Fig. 3. — Electron cooling time scales for our baseline model (a) as a function of electron 
energy right before the end of the injection period (same time as in panel 2d), and (b) as a 
function of time for a fixed electron energy of 7 = 10 3 . 
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Fig. 4. — Effect of a changing magnetic field on (a) the time-integrated vF v (fluence) spectra, 
(b) the X-ray light curves at 1 keV (multiplied by 100 for clarity) and 3.5 keV, and (c,d) 
hardness-intensity diagrams (HIDs). In the HIDs, the soft flux is the integrated 0.1-2 keV 
flux, the hard flux is the 2-10 keV flux, and the hardness ratio is the ratio of the two. The 
abrupt shape of some of the HID tracks is an artifact of the 5 approximations used for some 
of our spectral calculations. The respective magnetic fields are indicated in the legend; the 
other parameters are the baseline model values discussed in the text and listed in Table I. 
The vertical lines in panel (a) indicate the photon energies at which the light curves in panel 
(b) were extracted. 
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Fig. 5. — Effect of a changing luminosity of the companion star on (a) the time-integrated vF v 
(fluence) spectra, (b) the X-ray light curves at 3.5 keV, and (c,d) hardness-intensity diagrams 
(HIDs). The respective values of are indicated in the legend; the other parameters are 
the baseline model values discussed in the text and listed in Table 1. The vertical line in 
panel (a) indicates the photon energy of 3.5 keV for which the light curves in panel (b) are 
extracted. 
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Fig. 6. — Effect of a changing electron injection luminosity (see legend) on (a) the time- 
integrated vF v (fluence) spectra, (b) the X-ray light curves at 3.5 keV, and (c,d) hardness- 
intensity diagrams (HIDs). See Figs. 4 and 5 for more explanations. 
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Fig. 7. — Effect of a changing electron injection spectral index (see legend) on (a) the time- 
integrated vF v (fluence) spectra, (b) the X-ray light curves at 3.5 keV, and (c,d) hardness- 
intensity diagrams (HIDs). See Figs. 4 and 5 for more explanations. 
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Fig. 8. — Effect of a changing low-energy cutoff of the electron distribution, (see legend), 
on (a) the time-integrated vF v (fluence) spectra, (b) the X-ray light curves at 3.5 keV, and 
(c,d) hardness- intensity diagrams (HIDs). See Figs. 4 and 5 for more explanations. 
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Fig. 9. — Effect of a changing high-energy cutoff of the electron distribution, 72 (see legend) 
on (a) the time-integrated vF v (fluence) spectra, (b) the X-ray light curves at 3.5 keV, and 
(c,d) hardness- intensity diagrams (HIDs). See Figs. 4 and 5 for more explanations. 
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Fig. 10. — Effect of a changing viewing angle (see legend), and thus, of the Doppler boosting 
factor, on (a) the time-integrated vF v (fluence) spectra, (b) the X-ray light curves at 3.5 keV, 
and (c,d) hardness- intensity diagrams (HIDs). See Figs. 4 and 5 for more explanations. 
The values of the Doppler factor D corresponding to the various inclination angles are: 
D(5°) = 4.60, £>(20 o ) = 2.88, £>(35°) = 1.60, D(70°) = 0.583. 



